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Abstract 

pl,: 

Fokker-Planck equation with the velocity-dependent coefficients is considered 
for various isotropic systems on the basis of probability transition (PT) ap- 
i— i ' proach. This method provides the self-consistent and universal description of 

friction and diffusion for Brownian particles. Renormalization of the friction 
coefficient is shown to occur for two dimensional (2-D) and three dimensional 
(3-D) cases, due to the tensorial character of diffusion. The specific forms of 
^h *| PT are calculated for the Boltzmann-type of collisions and for the absorption- 

Q_i| type of collisions (the later are typical for dusty plasmas and some other sys- 

tems). Validity of the Einstein's relation for the Boltzmann-type collisions is 
analyzed for the velocity-dependent friction and diffusion coefficients. For the 
Boltzmann-type collisions in the region of very high grain velocity as well as it 
t*^- | is always for non-Boltzmann collisions, such as, e.g., absorption collisions, the 

Einstein relation is violated, although some other relations (determined by 
the structure of PT) can exist. The generalized friction force is investigated 
in dusty plasma in the framework of the PT approach. The relation between 
this force, negative collecting friction force and scattering and collecting drag 



o 

(N 
O 

O ! forces is established. 

The concept of probability transition is used to describe motion of active par- 
^ ^h * tides in an ambient medium. On basis of the physical arguments the PT for 

a simple model of the active particle is constructed and the coefficients of the 
relevant Fokker-Planck equation are found. The stationary solution of this 
, equation is typical for the simplest self-organized molecular machines. 

PACS number(s): 52.27.Lw, 52.20.Hv, 52.25.Fi, 82.70.-y 

a ■ 



I. INTRODUCTION 



Brownian dynamics nowadays is in the focus of interest due to the wide new fields 
of applications: physical-chemical systems, so-called active walkers, e.g., cells and other 
objects in biological systems, dusty plasmas with natural and artificial grains and many other 
systems. The characteristic property of such systems is the velocity dependent friction and 
diffusion coefficients. Existence of the Einstein relation and even the correct specific forms of 
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the Fokker-Planck equation for such systems are not still completely clarified. In particular, 
the attempts to use the Langevin equation as a stochastic basis for derivation of the Fokker- 
Planck equation lead to non-sign- valued result. The different forms of the Fokker-Planck 
equation, such as so-called Ito and Stratonovich [1-4] ones, appear. For the systems close to 
equilibrium Brownian particles keep stationary random motion under action of the stochastic 
forces, which are compensated by the particle friction and thus, the work produced by the 
Langevin sources is equal to the energy dissipated in course of the Brownian particle motion. 
This energy balance is described by the fluctuation-dissipation theorem in the form of the 
Einstein law. Obviously, the fluctuation-dissipation theorem and the Einstein relation can 
be violated in the case of non-equilibrium systems (even in the stationary case), in particular 
in the open systems. Starting from the classical Lord Rayleigh work [5] many studies of the 
non-equilibrium motion of Brownian particles with an additional (inner, or external) energy 
supply have been performed. In particular, such studies are of great importance for physical- 
chemical [6,7] and biological [8] systems, in which non-equilibrium Brownian particle motion 
is referred as the motion of active Brownian particles. The dynamical and energetic aspects 
of motion for the active Brownian particles have been described recently on the basis of the 
Langevin equation and the appropriate Fokker-Planck equation [9,10]. Possibility of negative 
friction (negative values of the friction coefficient) for Brownian particles was regarded, as a 
result of energy pumping. For some phenomenological dependence of the friction coefficient 
as a function of the grain's velocity the one-particle stationary non-Maxwellian distribution 
function was found. 

The traditional formulations of the non-equilibrium Brownian motion are based on some 
phenomenological expressions for the friction and diffusion coefficients. In particular, it 
means that deviations from the Einstein relation, as well as the velocity dependence of these 
coefficients are postulated and high level of uncertainty for application of such models to 
the real systems takes place. Recently we considered another situation, when the kinetic 
coefficients can be calculated explicitly on the basis of microscopically derived Fokker-Planck 
equation for dusty plasmas [11,12]. It was recently shown [13] that in the case of strong 
Coulomb interaction of highly charged grains in dusty plasmas, due to ion absorption by 
grains, the friction coefficient can become negative. The necessary criterion for negative 
friction due to ion absorption is found as T = e 2 Z g Zi/aTi > 1 (here Z g , Z { are the charge 
numbers for the grains and ions respectively, a is the grain radius, Tj is the ion tempera- 
ture). The appropriate threshold value of the grain charge was determined. The stationary 
solution of the Fokker-Planck equation with the velocity-dependable kinetic coefficients was 
obtained and the considerable deviation of such solution from the Maxwellian distribution 
was demonstrated. The physical reason for manifestation of negative friction in that case 
is clear: the cross-section for ion absorption by grain increases, when the relative velocity 
between the ion and grain decreases, due to the charge-dependent part of the cross-section. 
Therefore, for a moving highly- charged grain (r ^> 1) the momentum transfer from ions to 
the grain in the direction of grain velocity can be higher than in the opposite direction. 

In this paper we develop more general approach, based on probability transition, to 
simplify the Fokker-Planck equation and to calculate the velocity dependent friction and 
diffusion coefficients for the different systems. On that way we find the various forms of 
probability transition for the Boltzmann-type and for absorption collision integrals. The 
crucial peculiarity of the exact expressions for the mentioned coefficients follows from the 
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exact representation of these coefficients through the function of probability transition: 
-it is impossible to define the coefficients independently not only for the processes, which 
describe the systems close to thermodynamic equilibrium, when the Einstein relation is a 
priory valid, but also for the systems in which there is stationary, but non-equilibrium state, 
or for the systems far from equilibrium; 

-any rigorous approximate model of the Fokker-Planck equation have to be based on self- 
consistent expressions for the friction and diffusion coefficients, based on the PT. 

As an example we consider a wide class of open or far from equilibrium systems, where 
the Einstein relation is not applicable. For active particles the suggested consideration 
can be easily applied by construction of the probability transition on basis of the physical 
arguments. 



II. PROBABILITY TRANSITION AND VELOCITY-DEPENDABLE FRICTION 

AND DIFFUSION COEFFICIENTS 

The appropriate kinetic equation describing motion of Brownian particles in some 
medium with the momentum exchange may be written as 

^M = J p (P,f) = /dqMP + q,q)/ 9 (P + q,t) 

-w(P,q)f g (P,t)}, (1) 

where f g (P) is the distribution function of Brownian particles (grains) of the mass M. The 
elementary process is change of momentum of the grain P to (P — q). The probability 
transition w(P,q) in Eq. (1) describes the probability for grain with linear momentum P 
to lose the momentum q. The Eq. (1) has a form of master equation. In general the 
probability transition is the function of time itself. To simplify Eq. (1) for the processes 
with the momentum transfer q P we have expand the right side of Eq. (1) by q. The 
result of the expansion is the Fokker-Planck equation for grains: 



df g (P,t) _ d 



dt dp 



4(P)/«(P) + ^ (Sy(P)/ 9 (P)) 



(2) 



The coefficients -Aj(P) and ^(P), as easy to see by expansion of the Eq. (1), are 
expressed explicitly through the probability transition w(P, q) by the relations (e.g. [3,14]): 



A(P) = J dqft«;(P,q), (3) 
Bij(P) = ^J rfqg^(P,q). (4) 

Let us suggest that probability transition w(P,q) is a function of only two vectors P 
and q. It means, for example, there is no, let say, drift velocity of the media, surrounded the 
grain, as well as some inner vector inside the grain, which can influence on the probability 
transition. For that case the general structure of the coefficients Ai and Bij is evident: 

MP) = mP), B y (P) = ^fB l} (P) + {5 i3 - ^L) B ± (P), (5) 
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where (3(P), B\\(P) and B±(P) are the functions of modulus P. Let us consider at first the 
stationary case to understand the form of the Fokker-Planck equation and solution, when 
the friction and diffusion coefficients are the functions of grain's velocity. On this basis, 
in particular, the well known problem related to the Ito [1] and Stratonovich [2] forms of 
Langevin and Fokker-Planck equations [3, 4] can be solved and applicability of the Einstein 
relation for the various kinds of functions of probability transitions can be investigated. The 
results can be used also for unstationary case, if the initial distribution function is isotropic 
(what is not always valid, naturally). Then, taking into account isotropy of the distribution 
function / s (P,t), the Fokker-Planck equation is given by 



dfg(V) d 

dt dVi 
or in the equivalent form: 

df 9 (V) , „ d 



P*(V)V i f g (V) + ^(D ll (V)f g (V)) 



VdV 



dt 



= <* + 1 ^> 



nV)f 9 (V) + ^(D l} (V)f g (V)) 



VdV 



(6) 



(7) 



Here s is the dimension of the velocity space and the scalar functions of V are the same 
as ones of P, but expressed via the equality P = MV. We use above the velocity variable for 
grains V instead momentum P and the diffusion tensor D^iy) = M~ 2 Bij(V). We also use 
these notations below. The functions j3*(V) and D\\{V) are determined via the transition 
probability as 



POO 

D\\{V) 
D ± (V) 



V 2 



l - J d s q(PqMP,q), 
^ |rf s q(Pq) 2 W (P,q) 



p2 

2M 2 P 



2(s - l)M 2 P 2 



Jd s q [P 2 q 2 -(P q ) 2 }w(P, q ), 



(8) 
(9) 
(10) 

(11) 



where (Pq) is the scalar product in the velocity space with dimension s. Eq. (8) can be 
rewritten in the form: 



1 



2P 2 



dqw(P,q) s(P q ) 2 /P 2 -q 2 



(12) 



We see that the three scalar functions of P, determined by the different moments of 
probability transition, permit to find the coefficients in the Fokker-Planck equation. 

For the anisotropic velocity distribution function or for presence of the external fields, 
which do not change the friction and diffusion (what is not always valid, naturally), Eq. 
(2) can be rewritten as: 



dfg(V) 

dt 



_d_ 

dV. 



KP\V)f g (V) + ^-^(D {l (V)f g (V 



+ 5. 



d 



V 2 ) dV 



(D ± (V)f g (V)) 



(13) 
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For simplicity external fields are not included in Eq. (13). The useful equivalent represen- 
tation of Eq. (13) has a form: 



df g (V) _d_ 
dt dV 



Vi ( p . iv)+ ^) uv)+Dlt(v /y^ 



V 2 dV, 



(14) 



The stationary solution of the Fokker-Planck equation with the kinetic coefficients from 
Eq. (8,10) for the grain distribution function f g (V) is 



fg(V) 



c 



D\\(Y) 



cxp 



— / dvv 



Dn(v) 



(15) 



where C is a constant, providing normalization. As easy to see from the Eq. (6) for isotropic 
case the stationary (as well as non-stationary) Fokker-Planck equation with the velocity- 
dependent coefficients has well defined form and the question of " renormalized" friction 
coefficient is solved completely by the Eq. (8). The uncertainty in choice of the Fokker- 
Planck equation in the forms suggested, e.g., in [1,2] and [4], created by attempts to connect 
the Langevin and the respective Fokker-Planck equations by one-to-one correspondence, 
starting from the Langevin equation. The real structure of this renormalization, due to 
tensorial character of the diffusion £)y(V) as follows from the Eq. (8), permits to reformulate 
the problem: what must be the structure of Langevin equation for s-dimensional case to 
be relevant to the single-valued Fokker-Planck equation, based on the specific probability 
transition. Because we use, as the basis, transition probability we can establish validity 
or violation of the Einstein relation between the friction and diffusion coefficients directly, 
without usual suggestion of Maxwellian form of the static distribution function for Brownian 
particles, which is valid for the equilibrium state (when the Einstein relation is fulfilled a 
priory) . In particular, existence of the Einstein (or some different from one) relation between 
the momentum-dependable coefficients can be investigated. Correspondence between the 
Fokker-Planck and Langevin equations for s-dimensional case on the basis of PT approach 
will be considered in detail in a separate publication. Below we find the probability transition 
and investigate the various cases for the PT and Fokker-Planck equations. 



III. THE BOLTZMANN-TYPE COLLISIONS 

Let us consider the Boltzmann collision integral for two species of particles-light compo- 
nent (called below atoms) with the mass m and grains with the mass M ( m <C M) , which 
interact one with another (generalizations can be easily done). To find for such a process 
the PT function w s (P, q) it is enough, for example, to transform the part of the Boltzmann 
collision integral, describing the loss of grains in the phase volume dP near the point P, to 
the variables P and q, where q is the momentum transferred during the elementary act of 
collision between atom and grain. Than, comparing the result of transformation with the 
Eq. (1), for 3D case we find 
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Tfl f 

w s (P,q) = — do 

jjj J 



q 2 da f q 2 
2/j, | (qn) | do \2/j, | (qn) 



fn 



m p ml 
M ~J\2 



2 

q n 



v qnj 



(16) 



Here d o = sin x ^0 is the element of the space angle for scattering with the differen- 
tial cross-section da, f n is the distribution function of atoms, q = (p' — p) = (P — P'), 
q cos x — (q n )j vector n is the unit vector along the velocity V of atom after collision in the 
system of center mass for colliding particles and \x is the reduced mass. The values of the 
limiting angles for x above and below are usually equal Xmin = and Xmax = tt, except some 
special situations, when the integral over x diverges, as for example, for the purely Coulomb 
interaction, when the known cutting with < Xmin < Xmax < ?r is necessary. Taking into 
account the relation q 2 / (2 | (qn) |) =| (v — V) | Eq. (16) can be rewritten in the form 
more useful for applications: 



w s ( V, q) = J dCl J dv5[q + fiv — /x | v | n]x 

/„(v + V) |v| ^(|v|, X ). 
a o 



(17) 



We can use the representation (17) to obtain the useful general expressions for the friction 
and diffusion coefficients: 



MV : 



Jdft y*dv[ V - (| v | n- v) 



x 



da 



/»(v + V) |v|— (| v|,x)- 



(18) 



D||(V) = I^/ d "/ dv[V ' (|v| n " v )] 2x 

/n(v + V) |v| ^(|v| )X ). (19) 
a o 

The similar expression can be written for D±. When the distribution function of atoms f n 
has the Maxwellian form with the temperature T, the function u>(V, q) can be simplified, 
taking into account the inequality q< P. 

If | V | < i>t, where for the thermal velocity of neutral particles (atoms) we use the 
notation vt, the distribution function in PT can be expanded and we arrive, with respective 
accuracy of the order /x/M~m/M, at the expression for w s (V, q): 

Tfl ( Tfl \ ^/ 2 f f 

w*(V,q) = -— n n i ^-^\ J dQ, J rfv5[q + /iv-Ai|v| n]x 

/ mv 2 \ , . , T . da ., s 
eXP \^Tj 1 V| (v,V) do" (l V|,X) - (20) 

If atoms with the density n n are considered as the point particles and grains have the radius 
a it is easy to find from Eqs. (18-20) the values of the coefficients (3 and Df 
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P(V) = 8y^(m/M)a 2 n n v T , (21) 



\/2n 

D \\( V ) = 8^(m/M 2 )a 2 n n Tv T , (22) 

We see that the Einstein relation is fulfilled: 

D(V) — M' 1 T (3(V). (23) 

Calculations of (5 and D\\ were done independently on the basis of the appropriate PT 
function. 

Let us now consider under general condition q<C P the opposite case | V vt to solve 
the problem of validity of the Einstein relation for the velocity- dependable coefficients of 
the Fokker-Planck equation for arbitrary values of grain velocities [4]. In fact the answer 
can be found already for the particular case of atoms, scattering by grains, considering as 
the hard spheres. Even in this simple it will be shown, the Einstein relation is 

violated for high grain velocity. In the limit of extremely high grain velocity we can use the 
simplest approximation for the distribution function / n (v + V) = n n 5{y + V) and from Eqs. 
(11,17-19) obtain: 

(3{V) = (2m/3M)Tra 2 n n V, (24) 



D\\(V) = (m/M 2 )7ra 2 n n V 3 , (25) 
For D± by similar calculations we find 

D ± (V) = , (26) 

Therefore in that limit instead the Einstein relation we find another relation between 
velocity-dependable /3(V) and D\\(V): 

D{V) = ^mV 2 (3{V). (27) 

This situation can be classified as far from equilibrium. For the considering case of uncharged 
spherical grains the simplest interpolation relation between the friction and diffusion coeffi- 
cients can be suggested: 

*V)-£(l + *£)«V). <*> 

The similar interpolations as well as the exact velocity-dependent relations for the arbitrary 
cross-sections can be found from Eqs. (18), (19). As it follows from the relations considered 
above /3*(V) ~ /3(V) for all the cases with accuracy m/M. In general, by use the expansions 
to higher order of the small relation q/P, the expressions for the functions (3*(V), D»(V) 
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and D±(V) can be calculated and possibility to neglect the difference between /3*(V) and 
f3{V) can be established for the considered Boltzmann-type collisions. 

The essential influence of velocity dependence on the values of the friction and diffusion 
coefficients and violation of the Einstein relation for the Boltzmann type collisions takes 
place only for the extreme velocities much higher than the thermal velocity of the light 
particles. For the case of low grain velocity on basis of general representation (18), (19) for 
the friction and diffusion coefficients it can be shown applicability of the Einstein relation 
for arbitrary cross-section of scattering. 



IV. ABSORPTION COLLISIONS 

Now turn to other type of collisions, namely to the absorption collisions, which are typi- 
cal, for example, for dusty plasmas and some other open systems. As well known the process 
of grain charging by absorption of the electrons and ions by grains leads to the stationary 
(but non-equilibrium) state in plasma discharges. In so called OML (orbital motion limited) 
approximation the electrons and ions approaching to the grain on the distance less than 
the grain radius a are assumed to be absorbed. It is clear that the absorption collisions 
cannot be described by the Boltzmann-type collision integral. The appropriate correct form 
of the absorption collision integral have been postulated and applied in [14,15]. The rigorous 
kinetic theory of the electron and ion absorption in dusty plasmas, which exists in parallel 
with the usual processes of electron and ion scattering by grains, have been developed in 
[11,12], where also the Fokker-Planck equation for the charged grains was justified. Below 
we use the simplest form of the Fokker-Planck equation for grains with a fixed charge (dis- 
tribution by charge assumed narrow, what is often in reality). Our aim here is to find the 
probability transition function for absorption and to demonstrate efficacy of such approach. 
More complicated cases can be considered similarly. We also ignore increase of the grain 
mass [16, 17], assuming that neutral atoms generated in the course of the surface electron- 
ion recombination escape from the grain surface into a plasma. Naturally this process also 
changes the momentum balance for the particles, but we will not consider this process in 
our present model. Then the kinetic equation for grains can be written as: 

E [dpf a (p){W a (p,P-p,Q)f g (P-p,Q) - W a (p,P,Q)f g (P,Q)}, (29) 

where / Q (p) is the distribution functions for the electrons (a = e) and ions (a = i) and 
Q = eZ g . The elementary process is the absorption of electron or ion with the mass m a . 
The probability of absorption is given by: 



W a (p,P,Q) =a c 



p 


_P 




p 


P 


M ~ 


m 


M ~ 


m a 



(30) 



where a c (v,Q) with v = (| ^ — R |) is the cross-sections for absorption (or collection) of 



M 

the light plasma particles by grain in OML theory: 
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* C (Q, v) = ™ 2 (l - ^4-) 9 (l - ^4-) . (31) 
V m a v z a/ V m a v z a/ 

Below on the basis of the Eqs. (29)-(31) two different, but related problems are solved. At 
first, to consider the problem of absorption in simplest form, we accept the next simplification 
in spirit of [16,17], namely we will consider absorption of neutral atoms of a one sort. It 
means we put Q = and instead summation by a save only notation with the index n (for 
the neutral atoms) in the Eqs. (29)-(31). Generalization of this simplest model to the case 
of charge absorption with a fixed charge is quite simple. 

The second problem in focus of our interest is connected with a real dusty plasma, 
when there is ion stream and so called drag force, applied to the grains and created by ion 
absorption and ion scattering, exists. 

Let us start with a system of neutral particles, absorbing by grains. The momentum 
transferred to the grain due to absorption, is equal to the momentum of the atom p colliding 
with the grain. Therefore the probability transition w c (P, q) for the considered case can be 
immediately found by comparison of the Eqs. (1), (29) 



w c (P,q) = /n(-q)o- c 



M m 



M m 



(32) 



If we choose the Maxwellian distribution for atoms and suppose, to consider the simplest 
case, that absorption is purely geometrical a c = na 2 (the particular case of the Eq. (31) for 
Q — 0), we easily obtain with accuracy ~ m/M, that (3*(V) ~ /3(V). A simple calculation 
leads to the relation between f3(V) and D\\{V): 

D\\{V) = 2M- 1 T(3(V) (33) 

This relation is different from the Einstein one already for low V and for V = coincides 
with the result obtained in [12,17], as the limiting case (Q = 0) of the Fokker-Planck equation 
for dusty plasma in the case of the dominant absorption collisions. Here we found this 
relation, on basis of the theory with velocity-dependent coefficients, based on the probability 
transition approach, developed above. For more general form of the probability transition, 
which contains the free functions ip(P),x(q) and a small parameter ( <C 1, w(P,q) = 
if){P) 0(|q + CP|) x(q) it 1S possible to show, that as above (with accuracy to ~ m/M) 
f3*(V) ~ f3(V) and the relation between f3(V) and D\\(V) also has a form independent from 

V and different from the Einstein one 

(3(V) = 2C M 2 J d qqX (q) (ggg) 

D\\(Y) I dqq* X (q) 4>{q) ' ( j 

This consideration shows that the structure of the Fokker-Planck equation for the pro- 
cesses, based on the Boltzmann-type collision integrals, is very different from the processes 
of other type, when the Boltzmann-type collisions are not relevant. For the first type of pro- 
cesses the Einstein relation is valid, even for the case of velocity dependent the friction and 
diffusion coefficients, but in the limit of low grain velocity | V \<^ vt ■ For non-Boltzmann 
type of momentum transferring the fluctuation-dissipation theorem does not exist even for 

V = 0, through some relation between the friction and diffusion coefficients can take place 
(specific for the each type of PT). These results have deep consequences for many physical 
systems, as well as for the systems of biological nature, e.g., cells, moving in solutions and 
other, so called, active walkers. 
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V. FRICTION AND DRAG FORCE IN DUSTY PLASMA 



In this section as an application of the developed theory we consider the problem of 
momentum transfer from the ion stream to grains in dusty plasmas. Due to ion absorption 
and scattering by grains the drag force, acting on grains, appears. This force plays a crucial 
role in many experimental observations and, probably, is important for formation of voids in 
dusty plasmas for both "earthly" [18-20] and microgravity [21] conditions; the theory of voids 
essentially based on the drag force was given in [22,23]. The ion drag force Dj consists of 
two parts, so called collection D Ic and scattering D/ s ones. In the paper of Barnes et al.[24] 
the approximate analytical expressions for the drag collection and scattering forces were 
done. Later on the theory of drag force was actively developed analytically and numerically 
[23-28] to improve the description of drag force. It was achieved for the scattering part D/ s , 
in particular, in the recent publications [27,28]. 

Here we are focusing on generalized description of friction and relation between the 
friction and drag forces. Let us calculate the generalized expression for the friction force Ff 
on basis of the Fokker-Planck equation for the charged grains ((Eqs. (3), (29)). We calculate 
here Ff only for ions, because we are interested the ion part of friction. Generalization for 
many species of the light components (electrons, atoms) is simple. By integration of Eq. (2) 
on momentum we find for the time evolution of the average momentum P g : 

d{U f t 9i) = ~ J dP A(P,y)/ 9 (P, t), (35) 

The function Aj(P, y) is a generalization of Ai(P) of the Eq. (5) for the case of existence 
of some additional vector in the probability transition. In the case under consideration this 
vector y is the momentum, determined by the velocity of the ion stream y = Mu: 

ii(P,y) = J dq^(P,y,q), (36) 

To find w(P, y, q), when there is ion flow in plasma, it is quite natural to use the shifted 
Maxwellian distribution function of ions. 

At first we consider the scattering part Ff s of Ff. The probability transition w s (P, y, q) 



in this case has the evident property: 

w s (P,y,q) = w s (P-y,q). (37) 
Then we obtain the general expression for the friction force Ff s : 

F,. = -f dP (P - y) (3(\P - y|) f g (P, t), (38) 

l(P - y) = (p z y)2 /^q[q-(P-y)]^(P-y,q), (39) 

A(P-y) = (P-y)/3(|P-y|) = J dqqw s (P - y, q), (40) 



In the limit cases of the scattering friction itself Ff 0s = Fj s and scattering ion drag itself 
Dis = Fj s the Eq. (38) can be rewritten as: 
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F% = -fdPPP{\P\)f g {P,t), P»y 



(41) 



F< = %y/ 3(|y|), P«y (42) 

For the momentum distribution function of grains f g = 5(P — Po) (if the force is calculated 
for one particle) Eq. (38) takes a form: 

F/. = -(P -y)/?(|Po-y|), (43) 

and describes both the friction force itself and the drag force. In general, according to 
Eq. (38), (43), there is competition between friction and acceleration. For the limiting case 
of the friction force itself F/os for the grain with momentum G and immobile ions and the 
opposite case - ion drag itself Dj s with the ion velocity u = G/M there is natural relation: 

F f0s ee F> (G) = -F /s <(G) ee D /S (G), (44) 

This picture can be easily generalized for a few species of the light particles. From Eqs. (17), 
(43) we find for Ff s the representation: 

Ffs = Tv^uff / rfv/i(v " v + u)[(u " v) ■ v] |v| Mv) ' (45) 



rXmax H (J 

a tr {y) = dSl — (1-cosx), (46) 

JXmin US' 

Here we use the limits for the angles of integration, taking into account to provide conver- 
gence for the Coulomb cross-section in the case of ion-grain scattering. Eq. (45) coincides in 
the limit case V = with the well known general formulae for the transferring of momen- 
tum from light to heavy particle in the process of scattering, which can be justified from the 
simple physical arguments, as it was done, e.g., in [29]. For the opposite limit case u = 
it describes the friction force Ff 0s for grain. This equation is also applicable, naturally, for 
the short-range scattering potentials, when Xmin — and Xmax — tt- The specific result for 
Dj s can be written for the Coulomb cross-section in the form: 

D /s = 2A Mur 2 lnA. (47) 

where A = y ^(m i /M)a 2 n i VT i and = ^Ti/rrii. The parameter Y = e 2 Z g Zi/aTi and 
usually T ^> 1 for dusty plasmas. The structure of the generalized Landau logarithm In A 
for dusty plasma is very important and have been recently considered in [27,28]. For very 
strong interaction the problem of the correct form of In A as function of plasma parameters 
is still not completely solved. 

Let us consider now the generalized collecting friction Fj c and in particular the collecting 
drag force D/ c . The formula of the similar to Eq. (45) structure, but with V = and with 
the collecting non-transport cross-section, instead transport scattering cross- sect ion, was 
applied also for the collecting drag force in [23,27,28] and other papers on phenomenological 
basis. Our goal here to investigate the collecting drag force by use the PT Eq. (32) for 
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absorption and to find the relation between the friction force and the friction coefficient for 
the collecting process. Recently in [13] the friction coefficient /3(V) in dusty plasmas was 
calculated explicitly for arbitrary grain velocity and parameter T. It was found, that (3(V) 
can change sign ("negative friction") from positive to negative for some velocity domain if 
the parameter T > 1. Here we reproduce the result of [13] for the total friction coefficient 
for a grain f3{V) (for the particular but important case r\ = (mjt> 2 /2Tj) <C 1 and arbitrary 



The terms in Eq. (48) are proportional to the atom density n a and to In A describe friction 
respectively with atoms and with ions by scattering. These terms are always positive. 
Other terms in Eq. (48) describe "negative friction" due to ion absorption by grains and are 
negative in the considering limit case 77 <C 1 if T > 1. "Negative friction" exists for small 
77, if the Coulomb scattering is strongly suppressed, when the Coulomb logarithm In A is 
small [13,27] (some additional reasons for it reduction are discussed in [28]), that is typical 
for strong interaction in dusty plasmas. The level of ionization has to be enough high to 
provide negative value of the friction coefficient. As we know these conditions in present 
are not reached in the experimental set-ups. Opportunity for manifestation of "negative 
friction" in the experiments requires as we already mentioned the special conditions. From 
Eqs. (35), (36) and (48) it follows straightforward that the collection friction force itself F/ 0c 
for a moving grain can be written as: 



To find the generalized collecting friction force Ff c (and therefore also the drag force) we 
have find the function w c (P, y, q) for collection. The crucial fact is that the relation similar 
to Eqs. (37) or (49) for PT function is not correct for the absorption in the case when an 
ion flow exists (u ^ 0), due to the different structure of the PT functions for the scattering 
and collection processes. 

To describe the ion stream with the velocity u 7^ we use again the shifted Maxwellian 
distribution of ions. As it easy to see in this case the PT function w c (P,y, q) can be 
expressed via w c determined by Eq. (32): 



For the momentum distribution function of grains f g = S(P — P ) (if the force is as above 
calculated for one particle) the friction F f c takes the form: 



r): 




(48) 



F /0c = -P O /3(77,r). 



(49) 



w c (P, y, q) = w c (P - y, q + mu). 



(50) 



Ffc = -A(P , y) = - /dq(q-mu) w c (P - y, q) 



(51) 



This equality can be written in the equivalent form: 



F /c = muAfl P - y |) - (Po - mu)/3 2 (| P - y |), 



(52) 
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where the coefficients $ are related with the zero and first moments of the PT function: 

Pi = J rfqw c (P - y, q) (53) 

= - 1 , 2 / dq [q • (P - y)] w c (P - y, q) (54) 
(P -y) J 

Let us consider the simple and practically important case, when the both vectors P and y 
are directed parallel or antiparallel to the same unit vector 1: P = 7Tol and y = y^X. Then 
for the friction we arrive at the following expression: 

F fc = 1 {u [m/3i(| tt - y Q |) + M(5 2 {\ tt - y |)] - 7r /3 2 (| tt - y I)} , (55) 

Here and below we use velocities related with the momenta 7r = MVq and y = Muq. 
Eq. (55) can be represented in the equivalent and explicit form: 



Ff c = -1 jmwo J ' 



t t- / j [v • 1 + wo - Vb] (i*o - V ) r I 
-mVb y dv — _ v fi[\ + l(u - V )\va c (v) \ , (56) 



Let us consider the special cases of Eq. (56). 

a) Vq ^> Mo, wo — > 0. In this case we arrive to the friction force F/ 0c : 



* , )c = -mV / rfv [V °' ( T V /2 Vo)] / t [v - Vo] va c (v), (57) 



This expression coincides with the collecting ion friction force, which leads to the negative 
collecting friction coefficient for T > 1 [13] and to the respective relations (48), (49). 

b) u ^> Vb, Vb — > 0. In this case the generalized collecting friction F/ c describes the 
collecting ion drag force Dj c : 

/u ■ v 
dv — — fi[v - u] va c (v) (58) 
XL 

This equation coincides with the expression for the collecting drag force, which have been 
suggested in [12] and used also in [27,28]. 

c) u 7^ 0, Vq 7^ 0. Temperature of ions is low | (m — V ) |^> Vxi, the ion distribution 
function tends to delta function n;<5(v + u — V ). Then Ff c tends to the force directed 
along the ion stream, which we denote as D/ 0c : 

D/oc = mriiU \ (u - V ) \ a c (\ (u - V ) |), (59) 

As it follows from Eq. (59) ion wind in the considered limit always accelerates grains. If 
uo is parallel to Vb and they are close one to another, but both are higher than vn, the 
enhancement of the drag force D/ 0c occurs with decrease of the relative velocity | (u — Vq) \. 
It is the consequence of the OML collecting cross-section, which probably can be observable 
for very fast grains, discovered in some experiments [30,31], or in cryogenic discharges. 
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VI. ACTIVE PARTICLES 



During the last decade investigation of motion of the self-organized objects, as, e.g., 
cells, is in a focus of interest, due to numerous measurements and observations of their 
dynamical behavior [8,32]. Our goal in this work to show, that construction of the relevant 
probability transition on basis of the simple physical requirements permits to justify the 
relevant description of such systems. In particular, we show, that for the simplest structure 
of PT for motion of an active particle, the known (and experimentally verified) structure of 
the velocity distribution of grains (cells) [32], which are able to the directed motion near a 
fixed non-zero velocity, can be justified. The coefficients of this distribution are calculated. 

Let us formulate some general conditions to find the structure of the PT for active 
particles. We suppose that the linear momentum, transferred from a grain (cell) to the 
surrounding medium, is created by loss of the inner energy of this grain. Below we ignore the 
processes of energy supply, which can be included separately in more complicated schemes. 

At first we assume that the transferred to medium momentum q < P is distributed 
near some fixed value qo- The frequency of generation of the transferring momentum q 
will be denoted as /i(q). It can be approximated, for example, by the product of the 
functions, describing distributions on modulus | q |= q and on the space angles 9, (p between 
the vectors q and P, if there are no others, besides q and P, characteristic vectors for 
the system. In this simplest case we can put /i(q) = u fii(q) fi 2 (9 , ip) , where u is the q- 
independent frequency of momentum generation by a grain (cell). Below we suggest that 
the (^-dependence of \x<i is absent. The distribution fii(q) can be gaussian or, for the limit 
case of very narrow distribution of the transferring momentum, can be approximated by the 
delta-function Hi(q) = ^oS(q — qo)/ q%, where A is a dimensionless constant. For the function 
Ii2(9) we suppose, that the angle 9 between the direction of the transferred momentum and 
the momentum of the active grain is enclosed between the values 71 — 9q < 9 < ir, where 9q 
is some acute angle. Due to this amplification of the grains takes place. Consideration for 
the 2D and ID cases is evident. 

The "weight function" A(9) for the angle 9 can be included to describe the axis- 
symmetrical non-homogeneity of amplification for the different angles 9. We can also include 
two additional weight "functions": S(P) and T(e — sq). The first one describes dependence 
PT from the modulus P of the momentum P, the second one provides that momentum 
transferring is possible only if the inner energy of a grain (cell) e is bigger than some fixed 
minimal value of the inner energy e , let say T(e — e ) = $(e — e ), where $ is the step-like 
function. 

Under these assumptions the PT for the momentum transferring w £ , due to loss of the 
inner energy of an active particle, can be written in the case under consideration as: 



The values of the friction and diffusion coefficients follows from general Eqs. (8)-(ll): 



w £ (P, q) = u T(e - e )E(P)A{9) t x 1 (q)^(9 - vr + 9 )$(n - 9). 



(60) 



P £ (V) = -27rz/ T(e - s )^- J dqq^M J 



d9 sin 9 cos 9 A(9), 



(61) 



7T-0O 




(62) 



7T-0O 
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To find the total expressions for /3(P) and D(p) we have add to the values (61), (62) of 
(3 £ {P) and D £ (P) the parts of the friction and diffusion coefficients aroused due to collisions 
between the cell, moving with the momentum P, and the surrounding particles (atoms) of 
the solution. For the case of 3D elastic collisions and the hard sphere interaction these parts 
were calculated on the basis of Eqs. (21)-(22). For the velocities of cells essentially less, 
than the characteristic velocity of atoms we can ignore the velocity-dependent multipliers 
and consider the parts of the coefficients /3 e i(V) = /3q and D e i(V) = Do, connected with 
the elastic collisions, as the constants. These constants, as well as the initial velocity- 
dependent functions (for V < v Ta ) (see Eq.(23)) are connected by the Einstein relation 
D el {V) = M~ l TI3 el {V) ~ M~ 1 Tf3 . 

If we make the natural assumption, that S(P) is a constant, what means that the PT 
w £ is not dependent on the cell velocity (as the process, which is determined by the inner 
state of the grains), we find from the Eq. (61), that f3 e (P) ~ 1/P- Due to this specific 
dependence the Fokker-Planck equation for cells can be written as 

d Jfl = ±{\Vfi, - K M V) + [D, + DfJgl] , (63) 

where K £ > is a constant, determined by the equality (3 £ = — ^ and the Eq.(61), and 
D £ is determined by the Eq.(62). For (3q can be used the Eq. (21). Finally, the result is 
Po = 8V2T mn a S / (3M^/7rm), where S = na 2 is the area of the grain and n a is the density 
of the atoms. The stationary solution of the Eq.(63) is the Gaussian distribution: 




f g (y) = Cexpl- 7 £-lV--±\ , (64) 



Here C is the constant of normalization and D% = (D + D £ ). 

The velocity dependence of the distribution function f g (V) (64) coincides with one, which 
has been found in [32] on basis of the phenomenological assumption, concerning the structure 
of the friction coefficient in the Langevin equation. This type of the velocity distribution 
function in our consideration is the consequence of physically clear choice of the perturbation 
transition function. It can be generalized for more complicated and practically important 
cases, when there are one or more (additional to P) vectors, which determine the direction 
of q. It can be some inner vector - "driver", which can be orientated on the external 
(e.g., surface) gradients of density, or temperature, or concentration of some ingredient in 
the ambient medium. In that case, naturally, the equilibrium state is not effectively one- 
dimensional. Active particle (e.g., cell) can turn during motion. These problems will be 
considered separately. 



VII. CONCLUSIONS 

Here we use the simple and effective way for concretization of the Fokker-Planck equation 
on basis of self-consistent determination of the friction and diffusion coefficients. Both 
are determined as the functionals of probability transition. This function possesses a very 
different structure for the Boltzmann-type collisions and the other ones. We found PT 
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for the Boltzmann-type collisions and proved that velocity dependent friction and diffusion 
coefficients are connected by the Einstein relation for the velocities of grain less than thermal 
velocity of the small particles. At the same time there is crucial violation of the Einstein 
relation for the higher grain velocity. Therefore in general for the velocity-dependable friction 
and diffusion coefficients even for the Boltzmann-type collisions the applicability of the 
Einstein relation is limited by not very high (nevertheless practically most important) values 
of the grain velocity. The velocity dependence of these coefficients and renormalization of 
the friction coefficient in 2D and 3D cases, as a consequence of the tensorial structure 
of diffusion are found. Because the Fokker-Planck equation is single-valued also for the 
velocity-dependent coefficients the problem of connection between Langevin and Fokker- 
Planck equation has to be reformulated as a problem of the relevant (to the Fokker-Planck 
equation) Langevin equation. 

For the non-Boltzmann collisions, e.g. for the absorption collisions, the structure of PT 
follows from the structure of the collision integral, obtained earlier [11,14,15] and leads, 
in particular, to the different from the Einstein's relations between the coefficients in the 
relevant Fokker-Planck equation already for the region of low grain velocities. 

As the example of application of the PT method to the more complicated systems we 
considered the generalized friction force in dusty plasma. The scattering and collecting 
parts of this force are determined by the generalized friction coefficient, as a function of 
ion stream and grain velocities. The remarkable fact is that the sign of the collecting 
friction coefficient can be negative for some plasma parameters, as it was recently shown 
[13]. Of cause realization of negative total friction coefficient for grains for dusty plasmas in 
experiment requires the special conditions, because other mechanisms of friction are exist. 
The ion drag force in such approach, as well as the friction force itself, are the particular 
cases of this generalized friction. The ion scattering and collecting drag forces are found and 
calculated for the various particular cases. Some phenomenological expressions, which have 
been used for calculations before, are rigorously proved and generalized. 

We also constructed the PT for the active particles (e.g., grains or cells) in an ambient 
medium for some simple situation. On basis of physically clear assumptions we found that 
the part of the generalized friction coefficient, responsible for self-motion, can possess the 
peculiarity jj, where P is the momentum of a grain. Given appropriate also the usual friction 
mechanism the stationary solution of the relevant Fokker-Planck equation is gaussian with 
a peak around some non-zero velocity. Some generalizations of the obtained results for more 
complicated cases are suggested. 

ACKNOWLEDGMENT 

I am very indebted to Sergei Filimonow for his kind and permanent support that pro- 
vided possibility to write this paper. 

Author thanks E.A. Allahyarov, W. Ebeling, U. Erdmann, M. V. Fedorov, V.E. Fortov, 
A.M. Ignatov, L. Schimansky-Geier, P.P.J.M.Schram, I.M . Sokolov, A.G. Zagorodny for valu- 
able discussions of the various problems, reflected in this work. 

Author is grateful to G. Morfill for invitation in the Max Planck Institute for Extraterrestrial 
Physics in Garching and to him and S.A. Khrapak for useful discussion of drag force in dusty 
plasmas. 



16 



REFERENCES 



[1] K.Ito, Mem.Amer.Mathem.Soc. 4, 51 (1951). 
[2] R.L.Stratonovich, SIAM J.Control 4, 362 (1966). 

[3] N.G. van Kampen. Stochastic Processes in Physics and Chemistry, Edited North-Holland 
Publishung, 1992. 

[4] Yu.L.Klimontovich. Statistical Theory of Open Systems I, Kluwer Academic Publishers, 
Dordrecht, 1995. 

[5] J.W. Rayleigh. The Theory of Sound, vol. I, 2nd edition (Dover, New- York), 1945. 
[6] P-G. de Gennes, Rev. Mod. Phys. 57, 827 (1985). 

[7] A.S. Mikhailov, and D. Meinkohn, In: Stochastic Dynamics, edited by L. Schimansky- 
Geier, T. Poschel, vol. 484 of Lectures Notes in Physics, (Springer, Berlin), 1997, p. 334. 
[8] M. Scheinbein, and H. Gruler, Bull. Math. Biology 55, 585 (1993). 
[9] F. Schweitzer, W. Ebeling, and B. Tilch, Phys. Rev. Lett. 80, 5044 (1998). 
[10] U. Erdmann, W. Ebeling, L. Schimansky-Geier, and F. Schweitzer, Eur. Phys. J. B15, 
105 (2000). 

[11] A.G. Zagorodny, P.P.J.M. Schram, and S.A. Trigger, Phys. Rev. Lett. 84, 3594 (2000). 
[12] P.P.J.M. Schram, A.G. Sitenko, S.A. Trigger, and A.G. Zagorodny, Phys. Rev. E 63, 
016403 (2000). 

[13] S.A. Trigger, and A.G. Zagorodny, Brownian Motion of Grains and Negative 
Friction in Dusty Plasmas, Phys. Rev. Lett., submitted (September 2002), 
http : / / arXiv. org /abs/ physics /02 1 1 1 3 . 

[14] A.M.Ignatov, Plasma Phys. Rep. 24, 677 (1998). 

[15] S.A.Trigger, and P.P. J.M. Schram, J.Phys. D32, 234 (1999). 

[16] S.A. Trigger, Contrib. Plasma Phys. 41, 331 (2001). 

[17] A.M.Ignatov, S.A.Trigger, W.Ebeling, and P.P. J.M. Schram, Phys.Lett.A 293, 141 
(2002). 

[18] D. Samsonov, and J. Goree, Phys. Rev. E 59, 1047 (1999). 

[19] H. Rothermel, T. Hagl, G. E. Morfill, M.H.Thoma, and H.M. Thomas, Phys. Rev. Lett. 
89, 175001 (2002). 

[20] R. P. Dahiya, G.V. Paeva, W.W. Stoffels, E. Stoffels, G. M.W. Kroesen, K. Avinash, 

and A. Bhattacharjee, Phys. Rev. Lett. 89, 125001 (2002). 
[21] G. E. Morfill, H. M. Thomas, U. Konopka, H. Rothermel, M. Zuzic, A. Ivlev, and J. 

Goree, Phys. Rev. Lett. 83, 1598 (1999). 
[22] J. Goree, G. E. Morfill, V. N. Tsytovich, and S. V. Vladimirov, Phys. Rev. E 59, 7055 

(1999). 

[23] M. R. Akdim, and W. J. Goedheer, Phys.Rev E 65, 015401 (R) (2001). 

[24] M.S.Barnes et al, Phys. Rev.Lett. 68, 313 (1992). 

[25] M.D.Kilgore et al., J.Appl.Phys. 73, 7195 (1993). 

[26] J.E.Daugerty, and D.B.Graves, J.Appl.Phys. 78, 2279 (1993). 

[27] S.A.Khrapak, A.V.Ivlev, G.E.Morfill, and H.M.Thomas, Phys.Rev E66, 046414 (2002). 
[28] S.A. Trigger, G.M.W. Kroesen, P.P.J.M. Schram, E. Stoffels, and W.W. Stoffels. Ion 

drag and Brownian motion in dusty plasmas. In: 29 EPS Conference, abstract and 

report, Montreux, June 2002, to be published. 
[29] E.M.Lifshitz, and L.P.Pitaevski. Physical Kinetics. Edited "Science" Publishing, 

Moscow, 1979. 



17 



[30] A. Melzer, A. Homann, and A. Piel. Phys. Rev. E 53 2757 (1996). 

[31] H.M. Thomas, and G.E. Morfill. Nature 379 806 (1996). 

[32] M. Schienbein, K. Franke, and H. Gruler, Phys.Rev E. 49, 5462 (1994). 



18 



